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ABSTRACT 

Constraints on an exact quintessence scalar-field model with an exponential potential are de- 
rived from gravitational lens statistics. An exponential potential can account for data from 
both optical quasar surveys and radio selected sources. Based on the Cosmic Lens All-Sky 
Survey (CLASS) sample, lensing statistics provides, for the pressureless matter density pa- 
rameter, an estimate of Omo = 0.311q 'Jj. 
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1 INTRODUCTION 

Recently, astronomers have suggested that the universe is roughly 
geometrically flat and accelerating its expansion. Usual types of 
matter, i.e. baryons, photons and neutrinos, fail to close the universe 
and generate attraction, which leads to a decelerated expansion. To 
reconcile observations with predictions of general relativity, two 
new components have been added to the energy budget of the uni- 
verse: pressureless cold dark matter (CDM) and dark energy, also 
known as quintessence, with negative pressure. 

Whereas the simplest explanation for dark energy is a cos- 
mological constant, one of the more interesting physical re- 
alizations of the quintessence is a cosmic scalar field mini- 
mally coupled with the usual matter action l Peebles & Ratra 1988 
Caldwell et al. 1998 1. Such a field induces the repulsive gravita- 
tional force dynamically, explaing the accelerated expansion of 
our universe, fc-essence, namely a scalar field with non canon- 
ical kinetic terms l Armendariz-Picon et al. 1999 1, models based 
on branes and extra dimensions, such as the Cardassian sce- 
nario jZh u and Fujimoto 2 004), and t he generalized Chaplygin gas 
<Kamenshchik et al. 200l[lZhu~2 004 1 can also drive an accelerated 
expansion. 

A flat model of universe with a subcritical CDM energy 
density Omo, with two- thirds of the critical density existing in the 
form of dark energy and with negligible amount in standard matter, 
matches observations from galaxy clustering (Bahcall & Fan 1998 
Caiiberg et al. 1998 1, large-scale structure (Peacock et al. 2001 
Verde et al. 2002 1, cosmic microwave background radiation 
<de Bernardis et al. 20001 IBalbi et al. 20001 IJaffe et al. 200T1 

IPryfa^Uih^002^type la supernovae iRies s et al. 19981 
Perlmutter et al. 1999 1, age of the universe and Lya forest 
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tHarun-or-Rashid & Roos 20011 |Wang et al. 2000) . So, to explain 
the bulk of the evidence, we require a universe filled, nearly 
exclusively, with exotic matter. 

Gravitational lensing, based on the quite simple and well 
understood physics of light deflection in a weak gravita- 
tional field, provides useful tools for constraining cosmologi- 
cal parameters. Gravitational lensing statistics l Cha e et al. 20021 
Waga & Miceli 1998| |Cooray & Huterer 1999| |Wang et al. 2000| 
Zhu*20ToT1jaTrret^L*2003r . effects of large-scale structure growth 
in weak lensing surveys (Benabed & Bernardeau 2001 1, Ein- 
stein rings in galaxy-quasar systems I Futamase & Yoshid a~2001l 
Yamamoto & Futamase 2001 1, clusters of galaxies acting as 
lenses on background high redshift galaxies l Seren o"2002l 
|Sereno & L ongo 2004i|, offer very promising opportunities to 
probe quintessence. Gravitational lens time delay measure- 
ments can also produce precise estimates of cosmological pa- 
rameters l Schechter 2004 1, and the effect of quintessence in 
such observations has been investigated l Giovi & Amendol a 20011 
Lewis & Ibata 2002 1. Results from techniques based on gravita- 
tional lensing are complementary to other methods and can pro- 
vide restrictive limits on the dark energy contribution, sometimes 
in disagreement with the concordance value I Kochanek 1996a l. 

In this paper, gravitational lens statistics is used to study a 
class of quintessence models with an exponential potential. Such a 
potential can drive a scaling solution at late time, i.e. the scalar field 
mimics the scaling of the dominant energy component. Usually, as- 
sociated with an exponential potential, a scalar field is considered 
such that its fractional density, fi^, is practically constant during 
part of the matter-dominated era. The usual view is that, Q, v being 
small at the beginning of the matter-dominated era due to bounds 
from nucleosynthesis theory, the exponential potential cannot pro- 
duce a today accelerated expanding universe Bean et al. 200TY 
This conclusion is mainly based on considering the equation of 
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state of the dark energy as nearly constant. Recently, Rubano & 
Scudellaro ( 2002 1 showed how, without such an a priori assump- 
tion, the exponential potential can reveal very interesting features. 
This new position in favour of the exponential potential relies on 
a very strong argument, i.e. on general exact solutions of the field 
equations, allowed by a suitable choice of the exponent of the scalar 
potential. It has been shown that, since the equation of state is 
not constant, the tracker condition cannot be treated in the usual 
way. Despite of its simplicity, this class of potentials is able to fit 
data from type la supernovae l Pavlov et al. 2002 1, galaxy cluster- 
ing ( Rubano & Sereno 2002 1 and cosmic microwave background 
<Di Domenico et al. 20021 . 

The paper is as follows. In Section|2| the quintessence model 
is introduced. Section|3|contains the basics of gravitational lensing 
statistics. Properties of the data sample used in the statistical anal- 
ysis are illustrated in Sec.|4| Section|5|lists the constraints on the 
cosmological parameters. Section [6] is devoted to some final con- 
siderations. 



2 MODEL DESCRIPTION 

A spatially flat, homogeneous and isotropic universe filled with two 
non-interacting components, pressureless matter (dust) and a scalar 
field ip minimally coupled with gravity, can account for the present 
day observational data. Here, we investigate a quintessence model 
with a potential of the form ( Rubano & Scudellaro 2002 1 



V(cp) oc exp{-Ay5>. 



(1) 



where A 2 = !2irG/c 2 . The exponential potential has been 
widely studied, in relation with both quintessence and infla- 
tion jPeebles & Ratra 19881 |Capozziello et al. 1996") . The particu- 
lar choice of A allows for general exact integration of the cosmolog- 
ical equations. I present only what is needed for the present work. 
The time dependent Hubble parameter H is 



H = 



2(1 + 2?) 



(2) 



3«(1 + t 2 ) 

and the pressureless matter density (pm), in units of the critical 
density p cr = 3H 2 /8nG, is 



l + r 



(1 + 2t 2 ) 2 ' 



(3) 



where t = t /t s is a dimensionless time and t s is a time scale of the 
order of the age of the universe. 

The relation between the dimensionless time t and the redshift 
2 is given by 

(4) 



i 2 (l + i 2 ) 

where to is the present value of t. This very simple cosmological 
model has two free parameters, t s and to, or, equivalently, Ho, the 
present value of the Hubble parameter, and Qmo, the present value 
of the pressureless matter density. As can be seen from Eq. 
JImo depends only on to. 

The angular diameter distance between Zd and z s in a flat uni- 
verse is, in units of c/ Ho, 



d A (z d , z s 



Hp 
H(z) 



dz. 



(5) 



where, now, H denotes the Hubble parameter as a function of the 
redshift. The dimensionless angular diameter distance in Eq. J5J 



depends only on Q,yio- For cosmological distances in an inhomoge- 
neous quintessence cosmology, I refer to Sereno et al. 1 2001 2002 1 



3 LENSING STATISTICS 

To calculate the statistics of gravitational lenses, simple assump- 
tions are made iKochanek 1993 Kochanek 1996a I. The standard 
approach is based on observed number counts of galaxies, on some 
empirical relations between velocity dispersion and absolute mag- 
nitude, and, finally, on the simple singular isothermal sphere (SIS) 
model for lens galaxies. 

The differential probability of a beam of intersecting a galaxy 
with velocity dispersion between a and a + da, at redshift z d in 
the interval dzd is 



d t driG / ) ( ) CC ^ 
dzdda da dzd 



where 



(6) 



s cr is the cross section for lensing event, and —rQ is the 
differential number density population. For a conserved comoving 
number density of lenses, Ug(z) = no(l + z) 3 . 

As usual in lensing statistics, all lenses are associated with 
optically luminous galaxies. The luminosity function (LF) of 
galaxies can be modeled by a Schechter function of the form 
jSchechter 19761 
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where a is the faint-end slope and n, and L, are the character- 
istic number density and luminosity, respectively. The parameter 
q characterizes the faint part of the LF, which is still uncertain; a 
value ~ — 1 implies the existence of numerous faint galaxies acting 
as lenses. 

Since early-type and late-type populations contribute to the 
multiple imaging in different ways, type-specific LFs are required. 
As a conservative approach, I do not consider lensing by spiral 
galaxies, as their velocity dispersion is small in comparison to E/S0 
galaxies. The luminosity of a galaxy is correlated with its line-of- 
sight stellar velocity dispersion a via the empirical relation 

L 



(8) 



where a, is the characteristic velocity dispersion. The exponent 7, 
for early-type galaxies, can be fixed to the Faber-Jackson value, 
7 = 4 I Fa ber& Jackson 19761 . 

Early-type galaxies are well approximated as singular isother- 
mal spheres. As shown in Maox & Rix 1 1993 1 and Kochanek 
( 1996a 1, radial mass distribution, ellipticity and core radius of the 
lens galaxy are unimportant in altering the cosmological limits. As- 
suming a flat model of universe, a typical axial ratio of 0.5 in a 
mixed population of oblate and prolate spheroids would induce a 
shift of ~ 0.04 in the estimation of fi M o I Mi tchell et al. 2004 1. 
well below statistical uncertainties. Since departures from spheri- 
cal symmetry induce a relatively small effect on lens statistics and 
the distribution of mass ellipticities is highly uncertain, spherically 
symmetric models supply a viable approximation. The cross sec- 
tion of a SIS is 

— " '($'(&)' 

where Dd, Dd B and D s are the angular diameter distances between 
the observer and the deflector, the deflector and the source and the 
observer and the source, respectively. 
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When the deflector population is modeled as SISs and it is dis- 
tributed according to the Schechter function, the differential proba- 
bility for a source at z s to be multiple imaged with image separation 
A8 by a galaxies between z d and z d + dz d is 

d ' T «M?) 4 (e)V«>" 



dz d dA8 



Hp 

H(z d ) 



D d D d , 



x exp 



V D ds J 
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(10) 



Afl 
AS. ■ 



where A6, = 8tt \*f) and A6 ; 

The total optical depth for multiple imaging a compact source, 
in a flat universe, is 

r(zs) = ^{(l + z B )d A (0,z B )f, (11) 

where F, = 167T 3 n» (-^) 3 (^) 4 T[l + a + 4/7]. 

The configuration probability that a SIS in a flat universe 
forms multiple images of a background source with angular sep- 
aration A9 is 1 Kochanek 1993 1 

d 2 T 



p c (A8)dA8 = - 



dz d dA6 
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Two corrections to the optical depth must be included: 
magnification bias and selection function. Magnification bias 
accounts for tendency of gravitationally lensed sources to be pref- 
erentially included in flux-limited samples due to their increased 
apparent brightness ITurner 19901 |Fukugita & Turner 1991| 
Fukug ita et al. 19 92 Kocha nek 19931 . The bias factor for a source 
at redshift z B with apparent magnitude m is given by 



B(m,z,M ) 



(dNsY 
V dm J 



(14) 



dN s 
dm 



(m + 2.5 log M, z)P(M)dM, 



Mo being the minimum magnification of a multiply imaged source, 
with value M = 2;P(M)dM = 2M 2 M~ 3 dM is the probability 
that a multiple image-lensing event causes a total flux increase by a 
factor M ( Kochanek 1993 1. The function dN s /dm is the differen- 
tial source number count in magnitude bins dm. The magnification 
corrected probabilities are 



p(m, z s ) = r(zs)B(m, z s 



(15) 



Lens discovery rates are affected by the ability to resolve mul- 
tiple source images (Kochanek 19931. Since observations have fi- 
nite resolution and dynamic range, a selection function must be 
included to correct the statistics for observational limitations. For 
the SIS model, selection effects can be characterized by the max- 
imum magnitude difference that can be detected for two images 
separated by A8, Am(A8), which determines a minimum total 
magnification Mt = M (f + !)/(/ - 1), where 2.5 log/ = 



Am jKochanek 1 993 1. Then, the corrected lensing probability is 
<Kochanek i993llKochanek 1996a> 

', /x f B(m,z,M{(A6)) , A/1S , A/1 
p(m,z) = p(m,z) / >-l y>Ap c (Ae)dA8 



(16) 



The sum of the lensing probabilities p for the selected sample gives 
the expected number of lensed sources, 



B(m,z,M ) 
t(z) I B(m, z, Mt(A8))p c (A8)dAe 



(17) 



where p (mi, Zi) = Pj. 

The corrected image separation distribution function, p c , de- 
scribes the probability that a source is lensed with the observed 
image separation. For a single source at z, it is (Kochanek 1993 
IKochanek 1996al 

p(m,z) B(m,z,M f (A8)) 



p'(m,z) B(m,z,Mo) 
B(m,z,M{(A8)) 



p c (A8,m,z) = p c (A8) 

= pJA8)- p " V '"''~'"" 1V " (18) 

v 'J B{m,z,M t (A8))p c (A8)dA8 

I can now introduce the likelihood function (Kochanek 1993 
IChae et al. 2 002 1, 



c=H(i~p[)1[pl 



(19) 



where Nl is the number of multiple-imaged sources and Ny is the 
number of unlensed sources, pi is the suitable differential prob- 
ability accounting for the whole of the data available for each 
lens system, i.e. the lens redshift and/or the image separation 
iChae et al. 20021 IMitchell et~ 2004 1. I use a uniform distribu- 
tion for the prior on the cosmological parameter Quo, so that, apart 
from an overall normalization factor, the likelihood can be identi- 
fied with the posterior probability. 



4 DATA SAMPLE 

Reliable cosmological constraints from gravitational lens statistics 
need an unbiased sample. Indeed, a statistical sample must be com- 
plete and selected with well-defined criteria. The statistical prop- 
erties of gravitational lensing strongly depend on the properties of 
the distributions of sources in redshift and in luminosity in the ob- 
servational selection wave band. 

Analyses of statistical lensing have been based either on op- 
tically selected samples or radio- selected sources. The two meth- 
ods present both advantages and problems. Radio surveys are suc- 
cessful in selecting homogeneous and complete samples of sources. 
Since a radio source is selected purely on its radio flux and spec- 
tral index and the radio properties of a lens are unaffected by op- 
tical properties of the lens galaxy or extinction in the lens galaxy, 
radio surveys can avoid some systematic errors that can affect re- 
sults based on quasar lens surveys at optical wavelengths. The main 
shortcoming with radio-selected samples is the lack of detailed in- 
formation on the global properties of their redshift and luminosity 
distribution. On the other hand, optical samples require problem- 
atic matching of independent data sets (Kochanek 1996a 1 and suf- 
fer problems concerning extinction, but have a full descriptions of 
the properties of individual sources. In order to take advantages of 
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Figure 1. Normalized redshift distribution of the flat spectrum sources in 
Marlow et al. 2001. The kernel estimator of the distribution (full line) and 
a Gaussian distibution with (z s ) = 1.27 and dispersion 0.95 (dashed line) 
are overplotted. 



the properties of observations in the two wave-bands, in this paper 
the statistical analysis will be performed with both the optical- and 
radio-selected samples. 

A proper modeling of the distribution of the lensing galaxies 
is central in lensing statistics. Different points of view have been 
expressed on the best choice for the distribution in the deflector 
population I Mit chell et al. 20041 IChae et al. 2004> . However, re- 
sults from several recent galaxy surveys containing a large num- 
ber of galaxies appear to converge towards concordant estimates 
I Chae 2003 1. Chae 1 2003 1 used data from the Second Southern Sky 
Redshift Survey (SSRS2) to derive relative LFs. The early-type LF, 
expressed in the B photometric system, has: a — —1.00 ± 0.09, 



where h is Ho in units of 



n* = (6.4 ± 1.9)xlQ _ V Mpc~ 
100 km s _1 Mpc -1 . The early-type characteristic velocity disper- 
sion can be fixed at a, = 192±34 km s" 1 iChae 2003L In this 
work, the SSRS2 LF will be adopted. 



4.1 The Cosmic Lens All-Sky data sample 

Radio-selected galactic mass-scale gravitational lens survey can 
provide very large, homogeneous, unbiased statistical samples of 
sources. The Cosmic Lens All-Sky Survey (CLASS, Browne et 
al. <2003ll2003l ; Myers et al. <2003l ) is the largest radio selected 
galactic mass-scale gravitational lens search project to date. Out of 
the about 16,000 imaged sources, including the predecessor project 
Jodrell Bank Very Large Array Astrometric Survey, a subsample of 
8958 flat-spectrum radio sources down to a 5 GHz flux density of 
30 mjy, with 13 lenses constitutes a well-defined subset suitable 
for statistical analysis I Browne et al. 2003 l. I use the data listed 
in Table 1 of Chae 1 2003 1 updated with the spectroscopic obser- 
vations by McKean et al. 1 2004 1, which identified the lenses for 
0445+123 and 0631+519 as early-type galaxies at z d = 0.558 and 
0.620, respectively, while that for 0850+054 is a spiral-type galaxy 
at z d = 0.588. 

I limit my analysis to the early-type lens galaxies, as in the 
standard approach to lens statistics (Mitchell et al. 2004 1. Whereas 
the description of the late-type galaxy population is plagued by 
large uncertainties, they contribute no more than 10-20% of the to- 
tal lensing optical depth. I exclude 02 1 8+357 and 0850+054, whose 
lenses have been certainly detected as spiral galaxies. There are also 
arguments for discarding 1359+154 which presents a multiple lens 
system with three deflecting galaxies l Mitchell et al. 2004 1. So, I 
perform the statistical analysis considering two cases separately, 



the first one with 1 1 lenses and the second one with 10 cases. I do 
not use the measured image separations of 1359+154, 1608+656, 
and 21 14+022, whose splitting are due to multiple galaxy systems. 

A proper lens statistics analysis requires an accurate knowl- 
edge of the global properties of flat-spectrum radio sources. The 
final statistical lens sample is well described at 5 GHz by a power- 
law differential number-flux density, \dN/dS v \ oc S^ v , with 
r\ — 2.07±0.02 l Chae 2003 1. When the source counts can be mod- 
eled with a power law, the bias factor for a SIS lens population 
reads 



B(m,z,M { ) 



3-t] V/ + 1 



(20) 



where / is the maximum detected ratio of the flux densities of the 
brighter to the fainter images. The final CLASS statistical sample 
has been selected such that, for doubly imaged systems, the flux 
ratio is < 10 and it is independent of the angular separation. 

Redshift measurements are only available for a restricted 
CLASS subsample. I model the redshift distribution of the un- 
lensed sources with a kernel empirical estimator <Vio et al. 19941 
|Ryden 1996} . Given a sample of N measured source redshifts, 
{z St i}, the kernel estimator of the distribution is 



9 



(21) 



where K is the kernel function. Since the redshift is limited to non 
negative values, I use the kernel 



K 



rcf [, %s 



K Ga . 
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where Kgbm is a Gaussian kernel 
1 



K Gau (x) 



-x-72 



,(22) 



(23) 



Using the kernel in Eq. 122> . the Gaussian tail extending to negative 
values of z a is folded back into the positive axis. The kernel width 
g is well approximated by g ~ O^^A^ -0,2 , with A the smaller 
of the standard deviation of the sample and the interquartile range 
divided by 1.34. 

I consider the spectroscopic observations of a sample of 42 
flat-spectrum radio sources from CLASS in Marlow et al. <200U . 
The 27 sources with measured redshift have a mean redshift of 
1.27, with a standard deviation of ~ 1. In Fig.Q the histogram 
of the source redshift sample and the kernel estimator are shown 
together with a Gaussian distribution also used to model the distri- 
bution 1 Chae 2003 1. For the unmeasured lensed source redshifts, I 
set z s to the mean source redshift for the lensed sources with mea- 
sured redshift, (z s )i cnsc d = 1.93. 

Finally, according to the CLASS selection criteria, the com- 
pact radio-core images have separations gretar than A9 m i n = 
0.3 mas. The probabilities that enter the likelihood must be con- 
sidered as the probabilities of producing image systems with sepa- 
rations > A# m ; n . 



4.2 The quasar optical sample 

As discussed above, radio lens surveys have little information on 
the intrinsic redshift distribution of the sources. Since the optical 
depth to multiple image scales as this can lead to a serious 
problem in statistical analyses of lens samples. Kochanek 1 1996b 1 
showed the strong correlation between the mean redshift of the 
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Figure 2. Expected number of lenses as a function of Qmo f° r the CLASS 
sample. Short-dashed horizontal lines are for the conventional 68.3% confi- 
dence limit in the case 1 1 events (horizontal full line) are detected in a given 
observation. 



fainter sources and the cosmological model, so that uncertainties 
in the redshift distribution, or equivalently in the radio luminos- 
ity function, lead to systematic uncertainties in the derivation of 
cosmological limits. Furthermore, since the fraction of the sources 
identifiable as quasars steadily drops as the flux density of the 
sources decreases l Mufioz et al. 2003 1, any estimate of the radio 
luminosity function needs to divide the sources into two popu- 
lations, quasars and galaxies, contrary to earlier studies. These 
shortcomings will be addressed with systematic programs to es- 
timate the redshift distributions of the fainter flat-spectrum sources 
I Mufioz et al. 2003> . 

So, it still turns out to be very useful to perform statis- 
tical analyses based on optical lens surveys. To perform the 
analysis, I use data from six optical surveys for gravitationally 
lensed quasars: the HST Snapshot survey iMaoz et al . 1993i . the 
HRCam survey jCrampton, McClure & Fletcher 1992), the Yee 
survey «Yee et al. 1993T the NOT survey Tjaunsen et al. 1995> . 
the FKS survey tKochanek et al. 19951 and the ESO/Liege survey 
jSurdej et al. 1993) . A total of 893 high luminous optical quasars, 
plus five lenses, is considered iKochanek 1996a l. Matching differ- 
ent surveys requires a close examination of the selection effects and 
reliability of the individual data sample. I refer to Kochanek 119931 
1996a l for an extensive discussion. 

The differential quasar number count, dNQ/dm, can be mod- 
eled as I Kochanek 1996a) 



Figure 3. Normalized likelihood function from lensing statistics as a func- 
tion of Qmo for the CLASS sample in the case of 1 1 lenses. 



dN Q 
dm 



oc 10 



-a^m — rriQ (z)) 



+ 10 



-b(m-m (z))\ 



(24) 



where the bright-end slope is a — 1.07 ± 0.07 and the faint-end 
slope is b — 0.27 ± 0.07. The break magnitude mo(z) evolves as 



{m + (z — 1), z < 1, 

mo, 1 < z < 3, 

m - 0.7(2 - 3), 3<z, 



(25) 



with mo = 18.92 ± 0.16 at B magnitude. Since lens surveys and 
quasar catalogues usually use V magnitudes, I adopt an average 
B-V colour of 0.2 tBahcall et al. 1992llMaoz et al. 1993i . Finally , 
I use two selection functions, as suggested in Kochanek 119931 . 
one for space observations and another one for all the ground based 
surveys. 




0.4 0.6 
Q.M0 

Figure 4. Normalized likelihood function from lensing statistics as a func- 
tion of Qmo for the CLASS sample in the case of 10 lenses. 



5 RESULTS 

Let us now perform a statistical analysis to determine the param- 
eter JImo, which enter the model described in Section [2] We first 
consider the radio sample and, then, the optical one. 



5.1 Radio sample 

Let us first evaluate the predicted number of lenses, Eq. 1 1 7i . In 
Fig- El I pl°t the predicted number of lensed sources, til, in the 




Figure 5. Expected number of lenses as a function of Qmo f° r the opti- 
cal sample. Short-dashed horizontal lines indicate the conventional 68.3% 
confidence limit in the case of 5 events (horizontal full line). 
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Figure 6. Normalized likelihood function from lensing statistics as a func- 
tion of flyo for the optical sample. 



5.2 Optical sample 

Let us first evaluate the predicted number of lenses, Eq. 1 17i . In 
Fig.[3] I plot the predicted number of lensed quasars in the adopted 
sample, uq . In the case 5 events are detected in a given observation, 
the conventional 68.3% confidence limit is 2.480 < uq < 8.382 
<Gehrels 19861 . It is n Q = 5 for Q M o = 0.60, with ft M o ~ 0.29 
at 68.3% confidence limit. 

The expectation value and the standard deviation of the like- 
lihood distribution, see Fig. [6] agree with the mode and the central 
68.3% interval around the peak, respectively. It is Qmo = 0.6±0.2. 
For Qmo = 0.6, 5.1 lenses are predicted. 

No clear upper bound can be derived from the analysis of the 
optical sample. As can be seen from Figs. [5] even the Einstein- 
de Sitter model, where Q.mq = 1 and dark energy is not present, is 
compatible with the data. 



CLASS sample. The calculation of confidence limits is based on 
standard equations derived from Poisson statistics (Gehrels 19861. 
In the case 11 events are detected in a given observation, the 
conventional 68.3% confidence limit is 7.734 < n s < 15.42 
iGehrels 19861 . It is n L = 11 for Qmo = 0.41, with 0.23 ~ 
fiMO ~ 0.64 at 68.3% confidence limit. In the case of 10 lenses, 
the 68.3% confidence limit, 6.891 < n L < 14.27 jGehrels 19861 . 
corresponds to 0.27 ~ Omo ~ 0.73; Wl = 10 for Qmo = 0.47. 

As shown in Mitchell et al. 1 2004 1, the likelihood can be ap- 
proximated as 



C ~ exp 



N z {z a )p (z s )dz a 



(26) 



The likelihood is maximized for Qmo = 0.28_ ' 13 for 11 lenses, 
see Fig. [3] or 0.31^q for 10 lenses, see Fig. [4] Uncertain- 
ties denote the central 68.3% interval around the mode. Since the 
likelihood functions are slightly asymmetric, I also provide the 
expectation value and the standard deviation of the distributions 
P'Agostini 2004) . I get Q M o = 0.31±0.13 for 11 lenses, or 
0.35±0.14 for 10 lenses. 

The finite sample size induces an error in the estimated red- 
shift distribution. From the initial distribution estimated from the 
original TV = 27 source redshifts, I created 1000 data sets each 
containing TV points. Each data set is then used to create a new 
kernel estimator for the redshift distribution. The likelihood analy- 
sis is then repeated using the new 1000 redshift distributions. The 
resulting maximum likelihood estimates of Qmo present a disper- 
sion of 0.08. I have also verified that, when modeling the redshift 
distribution with a Gaussian distribution instead of the kernel esti- 
mator I Chae 2003 , Mitchell et al. 2004 1, conclusions are really un- 
affected. 

The main uncertainty in the estimation of cosmological pa- 
rameters comes from uncertainties in the assumed parameters of 
the luminosity function which describes the lens population. I cre- 
ated 1000 sets of galactic parameters n*, a» and <r*. The values of 
the parameters are extracted from normal distributions centred on 
the best estimates of each parameter and with standard deviation 
given by the associated uncertainty. The likelihood analysis is then 
repeated for each set of galactic parameters. The resulting distribu- 
tion of maximum likelihood estimates has a scatter of ~ 0.3, which 
give a similar uncertainty in the overall determination of Qmo ■ The 
main contribution comes from the uncertainty in the typical veloc- 
ity dispersion a* . 



6 DISCUSSION 

Exponential potentials in quintessence cosmologies have 
been recently re-considered and offer interesting perspectives 
(Rubano & Scudellaro 2002 1. As an example, an exponential 
potential can entail acceptable values for the CDM density at the 
nucleosynthesis epoch. 

A suitable choice of the exponent allows, through a transfor- 
mation of variables suggested by the Noether symmetry approach 
jCapozziello et al. 1996) , an exact integration of the Friedmann 
equations. This is the main attractive feature of the model we have 
been considering in this paper. 

Scenarios alternative to a cosmological constant have been al- 
ready tested with gravitational lensing statistics. Flat universes with 
dark energy with a constant equation of state has been considered 
(Waga & Miceli 1998|IChae et al. 20021 . Such an ad hoc modeliza- 
tion of the dark energy helps to discriminate between a cosmolog- 
ical constant and an evolving quintessence but it does not account 
for a general treatment of the dark energy component. On the other 
hand, physically motivated scalar fields deserve particular atten- 
tion. The model described in Section |2| presents, in my opinion, 
some very interesting properties. In general relativity, it is quite ex- 
ceptional to deal with exact general solutions. This enables a com- 
prehensive analysis and a full treatment of the model and of its 
phenomenological properties. 

Apart from this general argument in favour of exact solutions 
in general relativity, the examined exponential potential presents 
a further appeal from the statistical point of view. Since it has 
only two free parameters, i.e Qmo and Ho, following Bayesian 
arguments, it has the same a priori probability of a flat uni- 
verse with a cosmological constant and must be preferred to a 
flat quintessence cosmology with a constant, but yet undeter- 
mined, equation of state of the dark energy, wx ■ Despite of 
its simplicity, the examined exact solution reflects many proper- 
ties of the observed universe, with the same statistical confidence 
of a constant A-term <Di Domenico et al. 20 02 Pavlo v"et al. 2002( 
Rubano & Sereno 2002 1. It is compatible with the supernovae data, 
whose analysis gives Qmo = 0.15j_o'ois at 68% confidence level 
<Pavlov et al. 20021 . The peculiar velocity field and the perturba- 
tion growth of structure can also be accounted for in this scenario 
(Rub ano & Sereno 2002 1; data from galaxy clustering suggests a 



1 These models are characterized by three parameters: Ho, f^MO and Bx. 
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value of SImo = 0.18 ± 0.05 at 68% confidence level. The con- 
sidered model, not including radiation, is realistic only in the con- 
temporary or, at most, the recent past regimes of the life of the 
universe. However, a preliminary analysis of the anisotropies of 
the cosmic microwave background also gives concordant results 
<Di Domenico et al. 20021 . 

In this paper, I have considered complementary constraints 
from gravitational lensing statistics. This analysis is based on dif- 
ferent physics phenomena and on independent observations and can 
strengthen conclusions on a cosmological model. An analysis of the 
CLASS data sample provides a best-fit estimate, at the 68.3% con- 
fidence level, of Qmo = 0.3lig'i4 in the case of 11 lenses. When 
10 lenses are considered, the value of the best estimate increases of 
Af^MO — 0.04. In addition to radio-selected sources, optical sur- 
veys of quasars have been employed to constrain Qmo • This analy- 
sis slightly prefers higher values of the pressureless matter density, 
f2ivio = 0.6±0.2. The simple fact that the results from the two in- 
dependent samples are compatible at the 1-cr level is encouraging 
and we can consider constraints on cosmological parameters from 
gravitational lens statistics as quite reliable. Furthermore, estimates 
from lensing statistics are also compatible with results from both 
supernovae and galaxy clustering data. 

A number of systematic effects can plague lens statistics anal- 
yses (Chae 2003 1. Whereas most of the possible sources of errors 
are well controlled, the main uncertainty is connected with the 
adoption of the galaxy LF. As seen in Sec.|5l the main contribu- 
tion to the error in the estimate of SImo derives from the uncer- 
tainties in the early-type galaxy LF. Furthermore, there is no con- 
sensus among observationally derived results for the early-type LF 
jLoveday et al. 1992) IChae*2 003 1. Actually, a reliable determina- 
tion of the LF is needed to firmly draw conclusions on cosmolog- 
ical parameters. Nevertheless, as noted in Chae 1 2003 l, whenever 
the total galaxy LF is reliably determined, conclusions on cosmo- 
logical parameters are accurate if the partition of the total LF into 
the specific LF is careful. 
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